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We present ab initio studies for different Ni nanocontacts and show changes in the conductance 
of such constrictions due to atomic rearrangements in the contact. In particular we consider a Ni 
particle and show that the magnetoresistance can change from a few to 50 % and can even reverse 
sign as a function of the contact area formed between the particle and the leads. 
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I. INTRODUCTION 

Nanocontacts in the limit of a few atoms represent the 
ultimate miniaturization of electronic devices and earned 
a lot of attention during the last years. In particular, fer- 
romagnetic nanocontacts, their electronic structure and 
transport properties are of great interest and were inten- 
sively studied. 

The development of different experimental techniques 
allows to fabricate metallic nanocontacts. In princi- 
ple, the procedures can be classified into mechanical 
and chemical methods. The common mechanical tech- 
niques to produce contacts with a diameter in the order 
of a few nanometers are scanning tunneling microscopy 
(STM)^ and mechanically controllable break junctions 
(MCBJ)^ ^. It is known, that mechanical based methods 
are efficient and can be applied to a variety of materials. 
But they are not suitable for materials which tend to 
spontaneous breaking. Therefore electro-chemical meth- 
ods like controlled electro deposition^ ^ are preferable 
and allow the fabrication of contacts with a diameter 
of about 1 mm down to a few nanometers. The deposi- 
tion of the ions is regulated by an applied voltage. The 
disadvantage is, these methods are hard to control for 
nanocontacts in the limit of a few atoms^. 

Most experiments characterize the contact by mea- 
suring the conductance t hroug h the nanoconstriction. 
Monovalent metals like CrPJ^, A^, ArP^, and NsP^ 
show integer values of the conductance quantum in the 
order of a few = 2e^ /h down to l^o- The conduc- 
tance is very sensitive against tiny structural changes in 
the contact region which is refiected in the conductance 
histogram. 

During the last years, ferromagnetic contacts consist- 
ing of Fe^, Co^ and Ni^ ^ ^ ii 19 j^g^Q^j^g yery popular. 
Ferromagnetic contacts offer a further degree of freedom 
since the spin degeneracy is lifted. In particular, the rela- 
tive orientation of the lead magnetization can be changed 
by external magnetic fields. As a result a domain wall 
is geometrically pinned in the constriction which changes 
the transport properties and can cause giant magnetore- 
sistance (MR) expressed by the MR ratio 

MR = ~ X 100% . (1) 
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Here gp denotes the conductance for parallel alignment 
of the magnetic moments in the leads and gAP is the 
conductance for systems with anti-parallel lead magneti- 
zation. 

The goal of the experiments was to reach large dif- 
ferences between gp and gAP^ i-e. large MR ratios. 
Among al l the magnetic systems Ni was most of all 
analyzecPH^'^'^^1'^. Unfortunately, the results have been 
strongly contradictory also in comparison to theory^Q^^^. 

Some experimental results showed a half-integer 
conductance without an applied external magnetic 
field^ w hile oth er experiments did not support these 
observation d^l-*^^ ! -^ ^. Also in presence of magnetic fields 
the situation was confusing. Ono et al}^ reported a 
switching for the conductance from 2 e^//i to e^//i at low 
fields. Other authors claimed no changes by magnetic 
field^^. Furthermore, magnetostriction might be an ad- 
ditional source of uncertaint}®'^. 

Some papers reported about extr aordinary large MR 
ratios in Ni nanocontact^^^E^lI^n^, while others found 
maximum values of 40 

Theoretical investigation^^^^^ did not confirm neither 
the half-integer conductance nor huge MR ratios. 

The aim of this paper is to investigate structural 
changes in Ni nanocontacts and their effect on the con- 
ductance. Besides tensile and compressive strain we con- 
sider changes of the contact area between a nanoparticle 
and the leads which can be realized by a rotating parti- 
cle. For this purpose we simulate a Ni particle consisting 
of four atoms situated at the corner of a tetrahedron 
in different positions with respect to the leads. For all 
configurations we assume collinear magnetic order and 
investigate the parallel and anti-parallel lead magnetiza- 
tion. A treatment of the non-collinear magnetic order is, 
in principle, possible but numerically very demanding for 
complex contact geometrie^^^^*^. 



II. METHOD 

Our electronic structure calculations are based on den- 
sity functional theory (DFT)^^ and their extension to 
spin-dependent DFT^. We use the Korringa-Kohn- 
Rostoker formalism (KKR|^^ in the framework of the 
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local spin density approximation (LSDA)^^ with the pa- 
rameterization of the exchange-correlation potential in- 
troduced by Vosko, Wilk and Nusaiil^. 

For spherical potentials in the atomic sphere approxi- 
mation and in cell-centered coordinates the Green's func- 
tion writes as 

G(i?; r + R» , r + R»' ) = <5n,„' VI ^ i?2 (r < ) " (r > ) 

L 

L,L' 

(2) 

In this notation, the angular momenta (/,m) are repre- 
sented by the index L and n labels the atomic sites. The 
wavefunctions (r) and ^ (r) are the regular and ir- 
regular solutions of the single-site Kohn-Sham equations. 
The abbreviation r> (r<) denotes the vector with larger 

(smaller) absolute value of r and r^ In Eq. [ij G^^^f{E) 
are the structure constants as a function of energy, atomic 
sites, and angular momenta. 

The Green's function of the nanocontact is calculated 
in two steps. First, the Green's function of the host G 
consisting of two semi-infinite leads separated by a vac- 
uum spacer was calculated self-consistently. The struc- 
ture of the leads was assumed to be ideal fee with (001) 
surfaces. 

Second, the Green's function of the suspended 
nanocontact G was calculated by a Dyson equation: 

+ ^ G^'i: {E) At?: {E) (E) . 

The difference between the transition matrix tY{E) of 
the real system and the matrix f^^E) of the host system 
for the single scattering potential on site n is denoted by 
M^{E). 

Baranger and Stone^^ showed that the conductance g 
at zero temperature can be calculated via 

= — X 

I I dsi-G+(Ef;ri,r2)Vi^2G-(^F;r2,ri)-ds2, 

^ Si ^ S2 

(4) 

where G^{Ef; ri, r2) is the retarded (+) or advanced (-) 
Green's function at the Fermi level and is a double- 
sided derivative with respect to r^. The integrals in Eq.|4] 
are taken over the surfaces and ^2 located in the left 
and right lead indicated by the surface elements dsi and 
ds2 in the ideal leads, respectively. 

Following Ref. 36, the conductance can be written in 
terms of the structure constants G^'^, (E) and the KKR 
current matrix elements 

J2,L'iEF) = ^1/''^ RUEf;v)WRI,{Ef;v)\ (5) 



TABLE I: Irreducible representations of the group Cav 
and the corresponding angular momentum characters up to 

Imax ^ 3. 



irreducible representation 


angular momentum character 


Ai 




A2 


dx'^ —y^ 


A2, 


dxy 


As 


Pxi dxz 


As 


Pyi dyz 



as 

-—X 

EE E (JLMEF)-Jl.1,dEF))^ 

n,n' L,L' L" ^L'" 

{Jl',l'"{^f) - J2^n^Lf{EF))Gl'^^,{EF)G2)?^^n,{EF) . 

(6) 

The summation over n and is restricted to atomic 
spheres in the planes and 82- According to Biittikei'^, 
the total conductance g can be understood as a superpo- 
sition of individual transmission probabilities through 
individual eigenchannels m of spin a 

5=^Er™- (7) 

cr,m 

Scheer et al^ extracted the individual transmission 
probabilities for the eigenchannels from current-voltage- 
characteristics of atomic sized contacts between super- 
conducting leads. In particular, a projection of the con- 
ductance onto the atomic orbitals at the central atom 
in the contact is usecP^. Bagrets et al.^^ developed a 
method for the decomposition of the conductance into 
eigenchannels based on a KKR-basis set depending on 
site and angular momenta which is applied in this paper. 

In the general case of propagating Bloch functions the 
conductance channels are classified by the point group 
symmetry of the nanocontact. In case of C/^y symme- 
try the eigenchannels are characterized by the irreducible 
representations Ai,A2, A2/, and A5. A projection of the 
irreducible representations onto angular momenta up to 
I max < 3 is given in Tab. |l| 



III. CONTACT GEOMETRIES 

For our calculations, we consider three different kinds 
of Ni nanocontacts: a single atomic contact (Fig. [l^), 
a few atom constriction (Fig. [l]3), and a small Ni par- 
ticle bridging two semi-infinite leads (Fig. [l]^). For the 
single-atom contact and the few atom contact we focus 
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FIG. I: Sketch of the different nanocontacts suspended be- 
tween macroscopic leads. 



TABLE II: Distances between the plateaus (Pl, Pr) and 
the outermost atoms on the left or right hand side in (001)- 
direction of the Ni particle for the (2-2), (l-l) and (3-1) con- 
figuration. 



configuration 


distance d±^L (A) 


distance d±^R (A) 


(2-2) 


2.85 


2.85 


(3-1) 


3.23 


2.21 


(1-1) 


2.41 


2.41 



atom plateaus (Pl, Pr)- To simulate strain the distance 
d between the plateaus is changed gradually while the 
single atom (C) is kept fixed in the middle. 

To study the influence of symmetry changes in the con- 
tact region we embedded a square (C) with Ni atoms at 
the corners in between the two tips (Fig.[l]3). The corners 
of the square (C) and the tips (Tl, Tr) form a structure 
which is similar to fee in (001) orientation. The num- 
ber of Ni atoms on the central square (C) is gradually 
reduced from 4 to 1 resulting in the conflgurations GI- 
GS (see Fig.[l]D). 

Finally, we suspended a Ni particle between the 
plateaus (see Fig. The Ni particle consists of four 
atoms sitting at the corners of a regular tetrahedron. 
The edge length is equal to the next-nearest-neighbor 
distance of our fee Ni. The center of gravity of the Ni 
particle is flxed in the middle between the plateaus and 
the distance between the plateaus is flxed at d=7A8A. 
The tetrahedron is assumed to be a freely rotating par- 
ticle. We study three different conflgurations: a two-by- 
two (2-2), a one-by-one (1-1) and a three-by-one (3-1) 
conflguration. Because of the flxed center of gravity of 
the particle, the distances between the plateaus (Pl, Pr) 
and the outermost atoms of the particle are changing as a 
function of the particular orientation of the tetrahedron. 
For simplicity we show in Tab. [IT] the corresponding spac- 
ings in (OOl)-direction ((i±,L, between the left or 

right plateau and the outermost atoms on the left or right 
hand side of the particle, respectively. 



IV. RESULTS 
A. Single atom contact and strain 



on changes in the conductance caused by strain and sym- 
metry separately. For the Ni particle both effects are 
superimposed. 

All nanocontacts are sandwiched between two Ni leads 
in fee structure with an experimental lattice constant of 
do = 3.52 A. The surface normal is (001). The flrst two 
contacts are suspended between two tips consisting of flve 
atoms continuing the fee lattice. 

For the single-atom contact (Fig. [l^) the two tips are 
merged resulting in a single contact atom between four- 



To study the influence of strain we use the single-atom 
contact (Fig. [l^). The ideal fee lattice with the bulk 
lattice constant do is the "reference" system. 

Furthermore, we deflne a relative length variation of 
the plateau-plateau distance as Ad/do^ with do the bulk 
lattice constant and Ad = d — do^ the deviation of the 
plateau-plateau distance from the "reference" system. 
The relative length variation between the plateaus was 
changed in the range of —10% up to +10% in steps of 
5%. 

To characterize the single-atom contact magnetically 
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FIG. 2: (color online) Magnetic moments inside the constric- 
tion for the single-atom contact in ideal bulk structure for 
parallel (P) and anti-parallel (AP) lead magnetization in com- 
parison to the bulk magnetic moment. 



the magnetic moments for parallel and anti-parallel align- 
ment of the lead magnetization are shown in Fig. [2] For 
the parallel alignment (P) , the distribution of the mag- 
netic moments is symmetric with respect to the central 
atom in the single- atom contact. The magnetic moments 
in the planes {Sr-1) and Sl (Sr) are very similar 

to the bulk moment rriBuik = 0-62 jiB- The magnetic 
moment on top of the surface Pl {Pr) is increased with 
respect to the reduced coordination number. The mag- 
netic moment at the central atom is again bulk-like. 

For the anti-parallel configuration (AP) a domain wall 
is pinned in the constriction. As a consequence based 
on a collinear treatment of magnetism the magnetic mo- 
ment at the central atom is suppressed (Fig. [2| because 
of frustration. 

The total conductance values of the single-atom con- 
tact as a function of s train (Fig - |3|i) are in good agree- 
ment with experiment j^ l ^ l ^ * -^ -^^. Tensile strain reduces 
the conductance whereas compressive strain increases the 
conductance. This behavior is intuitively clear since the 
overlap of the atomic orbitals is reduced or increased, re- 
spectively. The overlap of the wavefunction is in a simple 
tight-binding-like picture a measure for the "hopping" of 
electrons. The change of the conductance for P and AP 
configuration with strain differs. An explanation of the 
different decay rates will be given below in terms of eigen- 
channels. 

The spin-resolved conductance for P and AP configu- 
ration is shown in Fig.jSfD and c, respectively. The current 
of majority electrons changes slightly with strain whereas 
the current of minority electrons strongly refiects strain 
in the P configuration. This fact will be discussed in 
terms of eigenchannels. 

The spin-resolved conductance of the AP configuration 
is degenerate with respect to spin because of symmetry 
{dti ~ "^^^ spin-resolved conductance is nearly the 
same as g^^. 
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FIG. 3: (color online) Total conductance for (a) parallel (gp) 
and anti-parallel (gAp) lead magnetizations of the single- atom 
contact, spin-resolved conductance for (b) parallel and (c) 
anti-parallel lead magnetizations, and (d) MR ratio as a func- 
tion of the relative length variation. 



The MR ratio (Eq. [T]) is plotted in Fig. [3]i. The MR 

ratio is positive in the range of —10 % to ~ +7.5 % rela- 
tive length variation. The sign reversal corresponds to a 
crossing of the conductance curves for the parallel and the 
anti-parallel configuration. For a relative length variation 
of 10 %, the conductance value for the anti-parallel con- 
figuration is higher than in the parallel case which causes 
a negative MR ratio. The maximum value of 14.7 % cor- 
responds to a relative length variation of —10%. The 
positive MR ratios are qualitatively in agreement with 
experiments^. 

The investigated nanocontact belongs to the point 
group C^v The results of the eigenchannel decompo- 
sition as a function of strain are given in Fig. [4j majority 
and minority channels in the P configuration in Fig. [4^ 
and b, channels of the AP configuration in Fig. 

Ni is a ferromagnetic material with the configura- 
tion 4s^ 3(fi for the conduction electrons, which leads to 
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FIG. 4: (color online) (^/e/t^ Symmetry resolved transmission probabilities for (a) majority and (b) minority states for the 
P configuration, and (c) transmission probabilities for the AP configuration for the single-atom contact as a function of the 
relative length variation, (right) Current density distribution in the plane S'l-S for a relative length variation of Ad/do = 0%. 



2 + 8 = 10 (partially) open channels for an ideal one- 
dimensional nanowire. The single atom contact is not a 
one- dimensional structure and therefore even more than 
10 conduction channels could exist. The following discus- 
sion is, however, restricted to channels with transmission 
probabilities > 0.10. 

In the P configuration, three majority channels are par- 
tially open (Fig. 4^). The main contribution corresponds 
to the irreducible representation Ai. If we would ana- 
lyze the channel compositioon, we would find 5, p^^ and 
dsz^-r^ orbitals at the single- atom bottleneck of the con- 
tact (see Tab. [l|. These particular orbitals are oriented 
in (OOl)-direction, i.e. the current direction. The other 
two channels are compatible with the representation A5 
and consist approximately of atomic Px (py) and the dxz 
{dyz) orbitals at single-atom position (C), or linear com- 
binations of them. These orbitals are oriented perpendic- 
ular to the current direction. Therefore, small changes in 
the distance d cause strong changes in the overlap of the 
wavefunctions. The reduced overlap induces a reduced 
conductance. This is the reason for the slightly enhanced 
decay of the channels T2^ and Tg^ in comparison with 



On the right hand side of Fig. [4| the current den- 
sity contributions for different eigenchannels in the plane 
6'l-2 are shown. From the symmetry of the current dis- 
tribution one can see, that the channel t|^ is dominated 
by s and Pz contributions, and the channel T2^ (T3^) is 
dominated by Px (py) orbitals. 



For the minority electrons (see Fig. [4p), six channels 
are partially open. Their corresponding irreducible rep- 
resentations are Ai, A5, and A2/. The comparison of the 
current density distributions shows the strong influence 
of d electrons which is already expected from the local 
density of states at the Fermi level (not shown). The 
strong d character at the bottleneck results in a strong 
distance dependence, which is illustrated in Fig. ^p. For 
all channels with representations A5 and A2/ the decay 
of conductance is slightly enhanced by increasing strain 
compared to channels with Ai representation. The num- 
ber of conduction channels for the minority electrons is 
larger than for the majority states. The decay rates of the 
channel contributions for minority states are very similar 
to the decay rates of majority electrons. But the number 
of conduction channels for minority states is larger than 
for majority states which is the origin of the stronger de- 
cay of the minority conductance in P configuration with 
increasing strain (Fig. [sJd). 

For the anti-parallel case, the analysis of conduction 
eigenchannels has been performed, too. The resulting 
transmission probabilities are plotted in Fig. |4j3. Here, 
four channels per spin with the corresponding represen- 
tations Ai, A5, and A2/ are formed. The channel T4^ 
with representation A2/ is marginally open. For higher 
relative length variations, this channel is closed. 

A comparison of the conduction channels shows that 
the anti-parallel alignment of lead magnetization changes 
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the composition of the channels and their transmission 
probabiUties drasticahy. The relatively strong reduction 
of the conductance in the parallel case is mainly forced by 
d-like contributions in the minority spin, while the ma- 
jority contributions change only slightly. The channels 
in the anti-parallel case show the same behavior like the 
majority states in the P configuration. In particular, the 
decay rates of the individual transmission probabilities, 
the number of open channels and their contributions to 
the transmission are very similar. 

Finally, the larger decay of conductance in the minor- 
ity spin compared with the majority contributions in the 
P configuration favors the crossing (see Fig. |4^, b and 
Fig. [s})). The spin-resolved conductances, their magni- 
tude as well as the number of conduction channels for 
the anti-parallel alignment of lead magnetization are very 
similar to the majority contributions in the P configura- 
tion (cf. Fig. and c). Because of degeneracy in the 
AP configuration no crossing of the spin-resolved con- 
ductances is observed. The crossing of the total conduc- 
tances for P and AP configuration (see Fig. 3^) and the 
corresponding change in the sign of the MR ratio (see 
Fig. [sji) is caused by an imbalance in the number of con- 
tributing channels to the spin-resolved conductances for 
both magnetic configurations. 



B. Small Ni clusters 

In a next step we focus on the infiuence on symmetry 
on the transmission and therefore we replace the single- 
atom contact by a few atom constriction. 

Details of the contact geometry have been discussed 
in Sec. 



TABLE III: Comparison of the spin resolved conductance 
values of the unconstrained P configuration of the sinlge-atom 
contact with the Gl configuration. 



Ill 



The symmetry of the contact is gradually 
reduced from C4V to lower symmetry by reducing the 
number of atoms in the central plane (see Fig. ^p). Ex- 
cept for G3 where in comparison to 02 one additional 
rotation by 180° and two refiections with respect the two 
atoms in C leaves the contact invariant. 

Before we discuss the conductance as a function of con- 
tact geometry we will concentrate on the magnetic struc- 
ture. In the case of parallel lead magnetization the mo- 
ments are ferromagnetically ordered in the ground state. 
The magnetization profile is symmetric with respect to 
the central plane C As a consequence of a reduced coor- 
dination numbers the magnetic moments increase. This 
behavior can be demonstrated by considering the mag- 
netic moments at different positions in the contact. The 
moments in the layers Sl-^ and Sr-1 (see Fig[T]3) are a 
slightly enhanced (m = 0.63 /i^) compared to the bulk 
value {rriBuik = 0.62 /i^). The magnetic moments in the 
surface layers {Sl and Sr) have a value of m = 0.64 fiB- 
A further reduction of the coordination number causes 
enhanced magnetic moments inside the contact (P^, T/,, 
C, Tr, and Pr) up to a maximum value of m = 0.80 /i^- 
The case of anti-parallel orientation of the lead mag- 
netization is more complicated. We are looking for the 
lowest energy configuration under the constraint of the 



conductance 


single-atom contact 


few atom contact 




2.10 


2.20 


9U [eVh] 


2.43 


0.97 



TABLE IV: Spin-resolved transmission probabilities for the 
dominating eigenchannels of the unconstrained P configura- 
tion of the single-atom contact with the Gl configuration. 



channels 


transmission for 




single atomic contact 


few atom contact 




0.88 


0.71 




0.54 


0.74 




0.38 


0.12 




0.51 


0.13 




0.50 


0.29 



anti-parallel orientation of the lead magnetization and 
the assumption of collinear magnetic order. As a result 
an abrupt domain wall is pinned between left and right 
electrode. Except geometry 04, the domain wall is lo- 
cated between the central atoms (C) and the left or right 
tip atom {Tl^ Tr). For symmetry reasons both positions 
are equivalent. The domain wall of 04 is located at the 
central plane. The moments of the two atoms on the 
central plane are suppressed. The moments in and near 
the surface are higher than the bulk values. The mag- 
netic moments on the central plane are in the order of 
m ^ O.TfiB- 

The conductance values and the corresponding MR ra- 
tios as a function of contact geometry are given in Tab.[V| 
Two aspects determine the general trend. First, the 
larger the contact area the higher is the conductance. 
Second, the lower the symmetry of the constriction the 
lower is the conductance. 

Before we discuss the effect of symmetry reduction on 
the conductance of the few-atom contact let us compare 
gp of the unstrained single-atom contact with the few- 
atom contact 01. Both contacts have the same symme- 
try C^V' The conductance of the single-atom contact is 
however higher than the conductance of 01. The spin- 
resolved conductance values are compared for both con- 



tacts in Tab. Ill The majority contribution g^^ is nearly 
unchanged whereas the minority contribution is drasti- 
cally reduced which is caused by blocking of eigenchan- 
nels. From the eigenchannel analysis in Sec. |IV A we 
know that the minority contributions are mainly d-\ike 
and the majority contributions are s-p-\ike. In Tab. IV 
the transmission probabilities of both configurations are 
compared. 

The dominating channels of the majority electrons 
have Ai and A5 character. As discussed before A5 is 
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TABLE V: Conductance values and MR ratios for Ni nanocontacts consisting of small clusters in dependence on the geometrical 
configuration in the contact region. 



geometric 




parallel configuration 




anti-parallel configurat 


ion 


MR 


configuration 




9U [eVh] 


gp [e^h] 


9n [eVh] 


m [eVh] 


QAP [e'^/h] 


[%] 


Gl 


2.20 


0.97 


3.18 


1.18 


0.78 


1.97 


61.38 


G2 


1.28 


0.98 


2.27 


1.45 


0.99 


2.44 


-7.25 


G3 


1.23 


0.91 


2.14 


1.04 


0.71 


1.75 


22.22 


G4 


0.79 


0.73 


1.52 


0.57 


0.57 


1.15 


32.28 


G5 


0.92 


0.69 


1.61 


0.57 


0.69 


1.26 


28.21 



twofold degenerate. The transmission probabilities of 
both configurations have the same order of magnitude. 

The channels of the minority electrons have also Ai, 
A5 and A2/ symmetry. The magnitude of the transmis- 
sion probabilities is however reduced by more than a fac- 
tor of two for the few- atom contact. The d contributions 
are blocked in the geometry Gl since they have to pass 
two single-atom contacts successively which acts like a 
filter for the strongly localized d states. 

With the help of the eigenchannel decomposition we 
can explain the changes of gp, in particular g^^, under 
the reduction of symmetry. Going from Gl to G2 by re- 
moving one atom the C4V symmetry is destroyed. The 
invariance under rotation around the z axis is lost. The 
twofold A5 does not any more exist. The corresponding 
channels are blocked and the conductance is reduced by 
1^0 (see Tab.[v|). Further reduction of the symmetry can 
be discussed analogously going from G2 to G4 or G5. 
In contrast, the conductance for G2 is higher than the 
conductance for G3, although the symmetry is increased 
going from G2 to G3. In this particular case the reduced 
contact area formed by the number of central atoms de- 
termines the conductance. 

The changes of g^ under symmetry reduction are not 
that pronounced since the channels are already only par- 
tially open. 

In the anti-parallel configuration the inversion symme- 
try with respect to the central plane is lost because of the 
magnetic order. As a result the qap values are usually 
smaller than the gp ones. Incoming majority electrons 
are scattered into minority states and vice versa. New 
eigenchannels are formed. 

A further reduction of symmetry in the central plane 
is not as important as in the parallel configuration since 
the pinned domain wall which is situated nearby the cen- 
tral plane acts already as a filter and requires a node of 
the eigenchannel. In some cases this fact can even cause 
gAP > gp (see G2). 

The individual conductance values are reduced for the 
few-atom contact with respect to the single-atom con- 
tact. The MR ratios are of the same order of magni- 
tude in agreement with experiment^^, can become even 
higher or reverse sign. Extraordinary MR ratios as 
measured^^ are however not obtained. 



C. Ni particle (tetrahedron) 

In this paragraph we discuss the results for the Ni par- 
ticle (see Fig. 

From some experiments, it is known, that structural 
changes can affect the transport properties strongljEEl^. 
Especially, for nanoconstrictions formed by controlled 
electro deposition of ions, the geometrical configuration 
is still uncleai^ and applied external fields, in principle, 
could cause atomic rearrangements inside the contact. 

We assume a Ni particle that, in principle, can rotate 
freely between the leads caused by an applied bias or by 
interaction of the particle magnetization with an exter- 
nal magnetic field. In our investigations, we assume that 
the particle rotates with respect to a fixed center of grav- 
ity in the middle of the constriction. For simplicity we 
consider three possible orientations of the Ni tetrahedron 
(see Fig. [l]^). 

First we will discuss the magnetic configuration of the 
Ni particle as a function of orientation with respect to 
the leads. The magnetic configuration is characterized 
by the averaged magnetic moments of all four atoms in 
the tetrahedron. 

The tetrahedra 2-2 and 1-1 are contacts, which have 
similar symmetry properties (see Fig. [l]^). This is re- 
flected in the electronic and magnetic properties, too. 
For instance, the distributions of the magnetic moments 
in the parallel configuration are symmetric with respect 
to the mirror plane. The averaged moments of the tetra- 
hedra reach values of 0.70 /i^ for the 2-2 and 0.94 /i^ for 
the 1-1 orientation. 

The domain wall in the anti-parallel configuration is 
located inside the tetrahedron for both orientations. The 
net magnetic moments are zero in both cases. 

The 3-1 contact is asymmetric. The three atoms on 
the left and the one atom on the right hand side have 
different distances to the plateaus (see Tab, [ill). 
ferent distances between the atoms in the tetrahedron 
and the atoms in the plateaus (Pl, Pr) are the reason 
for an asymmetric distribution of the magnetic moments. 
In the parallel configuration, the average magnetic mo- 
ment is 0.76 5 while in the anti-parallel configuration 
the net magnetic moment of the tetrahedron is 0.77 /i^- 
The domain wall is located between the left plateau {Pl) 
and the three atoms on the left hand side of the particle. 
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TABLE VI: Spin-resolved and total conductance for different orientations of the Ni particle in parallel and anti-parallel 
configuration. 



orientation of the 




parallel configuration 




anti-parallel configuration 


tetrahedron 




gu [eVh] gp [e^h] 


9n [eVh] 


m [eV^] 9AP [e^/h] 


2-2 


0.93 


1.33 2.25 


1.05 


1.05 2.10 


3-1 


0.90 


0.98 1.89 


0.81 


0.70 1.51 


1-1 


0.95 


1.25 2.19 


0.76 


0.76 1.51 



The calculated and spin-resolved conductances de- 
pending on the orientation of the tetrahedron are shown 
in Tab. I VII The calculated conductance values are of the 
same order of magnitude as for the constriction discussed 
above. A remarkable difference is that the conductance 
for minority states is larger than for the majority elec- 
trons in the parallel configuration which is in contradic- 
tion to the former cases of few atom contacts. However, 
this behavior of the spin-resolved conductances was ob- 
served for the single atom contacts with relative length 
variation Ad/ do < 7.5 % (cf. Fig.jsjo). 

The conductance values in parallel and anti-parallel 
configuration are the highest for the 2-2 orientation of 
the particle. These relatively high values are caused by 
the number of bridging atoms per lead. Each contact is 
contacted by two atoms respectively. This implies a large 
"effective" contact area in comparison with the other ori- 
entations, where one atom of the tetrahedron is contact- 
ing the plateau(s) and symbolizes a "bottleneck" (see top 
in Fig.^). 

The conductance in the anti-parallel case is nearly the 
same as for the parallel alignment of the lead magnetiza- 
tions and the spin-resolved values for the conductance are 
degenerated. The degeneracy is caused by the symmetric 
position of the domain wall. 

The conductance of the P configuration for the 2-2 
orientation is similar to the conductance for orientation 
1-1. The total conductance decreases from 2.2b e^/h to 
2.19e^//i and is influenced by the geometry. In orien- 
tation 2-2 and 1-1, the distances between the outermost 
atoms in the tetrahedron and the atoms in the plateau 
are nearly the same. The reduction of the number of 
bridging atoms from two to one can cause a blocking of 
the eigenchannels. 

In the anti-parallel configuration the domain wall is 
located in the tetrahedron for the 2-2 and 1-1 orienta- 
tion. As discussed above, the net moment of the particle 
is zero. The whole constriction has a mirror plane ge- 
ometrically and a point symmetry magnetically. There- 
fore the spin-resolved conductance is degenerate. But the 
conductance values are reduced compared to the parallel 
configurations. For the 2-2 orientation the conductance 
is higher than for the 1-1 orientation which is again a 
result of the reduced number of bridging atoms. 

By considering another orientation of the particle, one 
can end up with the 3-1 orientation. Now, the contact 
has no mirror plane. There are three atoms contacting 



TABLE VII: Calculated MR ratios in dependence on the ori- 
entation of the particle to the electrodes. The corresponding 
conductance values are listed in Tab.lVll 



orientation 


MR ratio 


initial (parallel) 


final (anti-parallel) 


[%] 


2-2 


2-2 


6.91 


2-2 


3-1 


49.42 


2-2 


1-1 


48.77 


3-1 


2-2 


-10.38 


3-1 


3-1 


25.25 


3-1 


1-1 


24.71 


1-1 


2-2 


4.09 


1-1 


3-1 


45.49 


1-1 


1-1 


44.85 



the left plateau {Pl) and only one atom contacts the 
right plateau (Pr). The conductance in the parallel and 
antiparallel case is determined mainly by the increased 
distance between the left plateau and the left hand side of 
the particle (see Tab.[ll[). In Sec. Ill we showed that ten- 
sile strain reduces the conductance which can be clearly 
seen by comparing the conductance values of the 2-2 and 
the 1-1 with the 3-1 orientation of the particle. In the AP 
configuration the conductance is slightly reduced com- 
pared to the parallel alignment of the lead magnetization 
due to the presence of a domain wall. 

Furthermore, the threefold symmetry of the particle 
with respect to the z axis does not fit the fourfold lead 
symmetry and causes blocking of the eigenchannels. This 
is reflected in the spin-resolved conductances, g^^ is only 
slightly reduced since the spherical symmetry of the s-p 
channel does not suffer from the symmetry misfit, g^^ 
with pronounced d character is, however, strongly per- 
turbed by the symmetry misfit and is reduced with re- 
spect to the other orientations of the particle. 

Assuming that the Ni particle can rotate freely trig- 
gered by an external magnetic field a variety of MR ratios 
can be obtained depending on initial and final orientation 
of the particle with respect to the leads. As mentioned 
above for simplicity we investigated three orientations of 
the particle which allows for exactly nine possible geo- 
metrical combinations for initial and final state of parti- 
cles orientation (see Tab. VII). The range of MR ratios 
spreads from +50 % to -10 %. The maximum values 
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are in agreement with experimental results^^ bu t does 
not support the huge values reported elsewhere^ ^ I4i i5 | i7| 
Interestingly the calculations predict a sign reversal of 
the MR ratio from positive to negative values depend- 
ing on the particle orientation. Recently it was reported 
about reproducible large MR ratios (+50 % to —20 %) at 
zero applied magnetic field in electromigrated permalloy 
nanocontacts (Ni8oFe2o)^^- Similar results were found for 
pure Ni nanocontacts. 

V. CONCLUSIONS 

We presented ab initio calculations of the ballistic con- 
ductance for different Ni nanocontacts. It was shown that 
compressive strain increases the conductance whereas 
tensile strain reduces the conductance of a single-atom 
contact. 

Furthermore, it was shown that the conductance of a 
few-atom contact strongly depends on the contact area as 
well as on the symmetry in the constriction. Reduction 



of contact area and of symmetry in the contact region 
with respect to the leads reduces the conductance. 

For all cases we discuss the magnetoresistance and 
obtained values of up to 50 % in agreement with 
experiment^ ^ The gigantic MR ratios^^ are not 
confirmed by our calculations. However, we predict a 
mechanism that can cause negative MR ratios. The 
mechanism is related to symmetry reduction in the con- 
tact region with respect to the leads and can be realized 
by a rotating Ni particle where the rotation is simulated 
by the external magnetic field. 
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